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Abstract 

We consider resonance tunneling of electromagnetic waves through an optical 
barrier formed by a stop-band between lower and upper polariton branches. 
We show that the tunneling through this kind of barrier is qualitatively differ- 
ent from tunneling through other optical barriers as well as from the quantum 
mechanic tunneling through a rectangular barrier. We find that the width of 
the resonance maximum of the transmission coefficient tends sharply to zero 
as the frequency approaches the lower boundary of the stop band. Resonance 
transmission peaks give rise to new photonic bands inside the polariton stop 
band in a periodic array of the barriers. 
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I. INTRODUCTION 



The effect of tunneling is well studied in the context of quantum mechanics (see, for ex- 
ample, book Recently, tunneling of electromagnetic waves has attracted interest owing 
to experiments with evanescent electromagnetic waves, [@-|5| which are a direct analog of 
wave functions of tunneling quantum particles. These experiments provide an opportunity 
to experimentally study the time evolution of tunneling wave packets - a subject of long- 
standing controversy (see, for example, review articles So far, three types of optical 
barriers have been considered in the context of the tunneling experiments. Historically, the 
first experiment with evanescent modes was carried out by Bose in 1927, where a prism with 
a beam incident at an angle larger then the angle of total internal reflection was used as a 
barrier |5[] . The same idea was used by Yeh [|TU| , who considered resonant tunneling of elec- 
tromagnetic waves in a superlattice composed of alternating layers with indexes of refraction 
such that the incident angle was greater then the angle of total internal reflection for one 
layer and smaller for the other one. The majority of time-of-flight theoretical considerations 
and experiments have been performed for an undersized wave guide ||,[2,11|, and photonic 
band gaps. [|3]|| All of the considered optical barriers were artificial structures. However, 
there exist natural optical barriers, which are well known, but were never considered in the 
context of light tunneling. We are referring to "restrahlen" or "stop" photonic bands, which 
exist in many dielectrics and semiconductors in the region of polariton resonances. It is 
well known that in the absence of attenuation, the transmission coefficient of light with a 
frequency within the stop band decreases exponentially with an increase in the width of a 
slab through which it propagates. Therefore, a simple slab of a dielectric with polariton 
resonances can serve as a "naturally made" photonic barrier. Properties of this polariton 
barrier, however, might be considerably different from properties of other kinds of barriers 
because of the different dispersion law of electromagnetic waves tunneling through the po- 
lariton barrier. It is important, therefore, to consider the tunneling properties of polariton 
barriers. 

In this paper we consider tunneling of electromagnetic radiation through a polariton bar- 
rier in a steady state mode. The point of interest in this situation is the resonance tunneling 
through a two-barrier system. Resonance tunneling of electromagnetic waves in a super- 
lattice structure was previously discussed in Ref. (pi. A tunneling effect in Ref. |TTJ arose 
when the angle of incidence for electromagnetic waves exceeded the angle of total internal 
reflection for one of the layers constituting the superlattice. The frequency dependence of 
the imaginary wave number in this case coincides with that of the square-barrier quantum 
tunneling problem. The situation considered in this paper is considerably different. The 
imaginary wave number in a polariton barrier layer depends upon the frequency in a peculiar 
way, demonstrating a singularity near the lower boundary, ut, of the polariton gap. This 
dependence does not have an analog in quantum mechanical systems. This leads to a new 
pattern of resonance tunneling states. 

Extending our system from two-barrier to a periodic array of alternating transparent and 
barrier layers, we obtain a new kind of photon band inside the initial forbidden band. These 
bands arise from the tunneling states and inherit many of their properties. In particular, 
we show that the bandwidths of these bands are determined by the widths of the respective 
resonances. The band structure of 2-D and 3-D photonic crystals made of materials with 
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frequency dependent dielectric permeability was numerically studied in Ref. |fl^-|14}|. The 
authors of Ref. 0] computed bands of two and three dimensional structures composed of 
metal rods (spheres in the 3-D case) with dielectric permeability of the form e = 1 — uj 2 /uj 2 , 



where u) p is the frequency of plasma excitations in the metal. The authors of Ref. |14| dealt 
with a photonic crystal built of a dielectric with e similar to that used in the present paper 
[see Eq.fll]) below]. One of the interesting effects observed in both studies was flattening of 
the photonic bands below the frequency u p , in the case considered in Ref. [T2|], or below 
ujt for Ref. [jH|] . We show analytically that this effect in the case of the dielectric photonic 
crystal is a direct consequence of the singular behavior of the dielectric function at ujt- We 
also discuss the influence of absorption upon tunneling properties of the polariton barriers. 



II. RESONANCE TUNNELING IN A THREE-LAYER SYSTEM 

We consider propagation of electromagnetic radiation of frequency w in a system that 
consists of two parallel absorbing dielectric layers, each of thickness a, positioned along the 
x-axis and separated by a distance b. The dielectric permeability of the dielectric material 
is given by 

(u + il? ~ uj 2 l 

e{u;) - £oo- — ——2 -, {!) 

[oj + 27) — OJf, 

where 7 is the attenuation, and is the background dielectric permeability. Both layers 
are considered optically isotropic, so that TE and TM waves are not coupled. This allows 
us to neglect the vector nature of the electric field. We consider the case of normal incidence 
so that the electric field in our system can be presented as follows: 

E(x) = e ik ° x + re~ ik ° x > x > —00, 

= aie ikx + bie~ tkx a > x > 0, 

= a 2 e iko{x - a) + b 2 e- iko{x ~ a) a + b>x> a, (2) 

= a 3 e ik{x ' a ' b) + 5 3e -ifc(x-a-fc) 2a + b>x>a + b, 
= te ik °( x ~ 2a - b) x>2a + b, 

where k — uj/c is the wave- vector in vacuum, c is the speed of light in vacuum, and r, a 1; 
bi, a2, &2, CL3, bs and t are the complex amplitudes of the plane waves in each of the five 
different regions. The wave number k in the dielectric layers is determined by the expression 
following from Eq. ([!]): 



[uj + ry) 2 - lo\ 
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k = kri\ £00 ^ ' 7T5 "T- ( 3 ) 



Taking into account regular boundary conditions at each of the boundaries between dif- 
ferent materials, we obtain the following expression for the complex transmission coefficient 
t: 
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In the case of a lossless dielectric (7 = 0) the expression for the transmissivity of the system, 
T —\t 2 \, can be obtained as 



1 + + ^) 2 sinh 2 (/ta) cos 2 [K(a + b)] 



(5) 



where K(a + b) is given by 



cos [K(a + b)} = cosh(/ta) cos(A;o&) + 



1 / K 

2 [k~ 



sinh(/€a) sin(A;o&). 



(6) 



We assume here that the frequency of the wave falls into the region ujt < oj < ujl so that the 
wave number k = in is imaginary, and propagation of the wave through the dielectric layers 
is evanescent. Eq. (|5[) describes optical tunneling through a forbidden "band." However, for 
a set of frequencies such that 



cos [K(a + b)] = 0, 



(7) 



the transmission coefficient, T, is equal to 1, and tunneling through the system becomes 
resonant. Though Eq. (|6|) is similar to the respective equations of Ref. |I(J and of the 
quantum mechanics tunneling problem, the different frequency dependence of the wave 
number k causes different behavior of the solutions to this equation. 

In order to determine the number of resonance frequencies, it is convenient to rewrite 
Eq. (0) as 



1 fkn 



-) 

K k J 



tanh(/«a) = cot(k b). 



If boundary frequencies ujl and ojt satisfy the inequality 
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LU L - LO T < — , 




(8) 



(9) 



the respective interval of wave numbers form ujtI 'c to ujl/c can only accommodate less than 
one period of cot(A;o&)- In this case Eq. (|8|) has one solution if 

k atan(k b) > 2. (10) 

If inequality (^) does not hold, the number of resonance frequencies is equal to 

(uj l - u T )b 



N 



+ 1 
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if inequality (pi]) is satisfied, and it is iV — 1 otherwise. The transmission, T, as a function 
of frequency is presented in Figs.l for two different values of the widths of the layers a and 
b. One can see that as the ratio of a to b increases, so does the number of resonant peaks. 

Both plots in Fig. 1 show a similar pattern: the separation between resonant peaks 
remains approximately constant while the half-width T of peaks sharply decreases as the 
frequency approaches ujt- The resonant half-width, T, can be found by expanding the 
denominator of Eq. (|5|) near the resonance frequency u r up to the term quadratic in u — u r . 
If a resonance frequency, u r , is close to the lower boundary of the gap, Ut, one can find for 

r 



2^2 (5m) 2 
k Q a 



exp 



(12) 



where 5u = u T — u r , and = e{oj 2 L — oj t )/{2uj t ). Eq. (|T2"D describes an extremely sharp 
non-analytical behavior of the width of the resonance near the low-frequency boundary of 
the gap. 

To take absorption into account, one can replace the parameter 5uj by the following 
expression 



5u abs = (uj t ~ UJ r ) 2 + 7 2 . 



If cut — oj r < 7 the width of the resonance becomes proportional to exp (— co^/7) 1 / 2 and 
remains small for 7 < uo*. However, the width grows exponentially fast with an increase 
of the absorption rate, and the height of the peak decreases correspondingly. Therefore, 
the resonances nearest to lot are the most vulnerable with respect to absorption. One can 
see this in Fig. 2, which shows that these resonances disappear first when the absorption 
becomes larger. It also follows from this analysis that the width of the resonance peaks 
close to ujt is determined solely by absorption. This behavior of the resonances is specific to 



our particular model and does not take place in the case of Ref . |10[ , where the wave- vector 
is finite for any uj in the forbidden band. 

When u r approaches ujl, the wave-vector n(uj r ) becomes smaller and T increases with 
frequency, reaching at ojq ~ ujl a certain value, which depends upon all the parameters 
of the system. For these frequencies, absorption, when small enough, does not contribute 
significantly to the width of the maxima. 



III. PHOTONIC BANDS IN THE POLARITON FORBIDDEN GAP 

In the case of a periodic array of the three-layer sandwiches considered in the previous 
section of the paper, one has a set of propagating bands instead of single resonance frequen- 
cies. In this situation, Eq. (|) presents the dispersion equation of the propagating modes 
with K being the Bloch wave number. The resonant tunneling frequencies discussed in the 
previous section are the solutions of the dispersion equation at the center of the Brillouin 
band, K(a + b) = tt/2, where K changes between and tc within the band. It is clear, 
therefore, that each of the resonances gives rise to a corresponding band. This conclusion 
is supported by the dispersion curves plotted in Figs. 3. One can see from this figure that 
for each resonant frequency u) rn , where n is the number of a given mode, there is a distinct 
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branch of the dispersion curve. It is interesting to note that branches corresponding to the 
resonant frequencies close to lot, have much smaller dispersion than branches which corre- 
spond to lo t ~ u L . One can also notice that regular bands, which appear outside of the 
polariton gap, also become less dispersive for frequencies near lot- This observation agrees 
with the results reported in Ref. for two-dimensional systems. The advantage of the 

one-dimensional model is that we can analytically determine the cause for this behavior. 
Assuming that the width of a band, centered at u rn , is much smaller then the frequency 
itself, one can expand the right-hand side of the dispersion equation (|6]) with respect to 
lo — uj rn and obtain the approximate dispersion equation as follows: 

( K k \ 

lo - co rn = (-l) n r (uon) - — I sinh [K(uo 0n )a] cos [K(a + b)} . (13) 

\k k J 

If one combines this equation and Eq. flT2| ) for the resonance width T, the expression for the 
band width, Au n , for bands near the lower edge lot becomes 

A ^ 

LO* 

This expression also holds for frequencies at the pass-band side of lot and explains, therefore, 
the flattening of the photon dispersion curves near lot observed in Ref. [|I2|-|I4||. Both the 
band width Aoo n and the width of the resonance F tend to a certain finite value when lo 
approaches lot- The prefactor (— l) n in Eq. (Jf^) shows that the branches have an alternating 
sign of the dispersion: positive for the "even" branches and negative for the "odd" ones. 
This feature is clearly seen in Fig. 3. 



IV. CONCLUSION 



We consider a polariton stop-band (restrahlen region) as an optical barrier for tunneling 
of electromagnetic waves. We show that the peculiar frequency dependence of the tunneling 
penetration length of electromagnetic waves results in tunneling properties that are qual- 
itatively different from those of other optical barriers. We carry out a detailed study of 
the resonant tunneling of electromagnetic waves through a three-layer sandwich, where a 
dispersionless (vacuum) layer was placed between two identical layers with dielectric per- 
meability allowing for the polariton stop-band. It was found that the number of resonance 
frequencies depends upon the frequency width of the stop-band, lot — ^t, and the spacial 
width, b, of the vacuum layer. This number increases with an increase of lot — lot and b. 
The latter dependence appears paradoxical, since it implies that no matter how far away 
tunneling layers are from each other, they are capable of providing conditions for resonant 
tunneling. The situation becomes clear if one recollects that we deal with the steady state 
situation, where the flux of energy in the system is fixed by an external source. For a large 
system, a wave has to travel a long distance between the layers before it reaches the steady 
state condition. This fact is crucial for an experimental observation of resonance tunneling 
since one has to have a steady enough source of light and be able to maintain a complete co- 
herency of the wave while it travels between the layers. Let us consider, for example, GaAs, 
for which lot = 5.1 x 10 13 s _1 , and lot = 5.5 x 10 13 s _1 [|15|| . The vacuum wavelength in this 
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case changes between 34/im and 37//m when frequency sweeps over the polariton gap. If 
one makes the distance b between the layers equal to, say, 350/xm (approximately ten wave- 
lengths), then the number of observable resonances would be equal to 1 or 2 depending upon 
the width of the dielectric layers. In order to observe greater number of the resonance one 
could use materials with wider polariton gaps, for example, MgO, where ujt = 7.5 x lO 13 ^" 1 
and Ul = 14 x 10 13 s _1 In this case even for b as small as 60/zm, four resonance peaks 
could be observed. The actual number of resonance maxima also depends upon relaxation 
characteristics of the layers. Though these characteristics can vary even for different sam- 
ples of the same type of material, one can consider 7 = 0.02 [T6] as a typical value and 
use it for quantitative estimates. It is seen from Fig. 2 that the relaxation of this order of 
magnitude could reduce the number of observable resonances by at most one maximum if 
the configuration of the system allows for a maximum near lot- More specific conclusions 
regarding the number of the expected resonances could be drawn from the results of the 
paper for each particular experimental configuration. 

The width of the resonance peaks of the transmission coefficient was found to have a 
sharp nonanalytical frequency dependence at the lower boundary of the polariton gap, ojt 
and to saturate at a certain finite value when approaching the upper boundary, ujl. This 
fact results in different reactions to absorption for resonances occurring in vicinities of ujt 
and ujl . In the first situation even a small absorption completely determines the real width 
of the resonance, and these resonances are washed out first with an increase of absorption. 
The resonances in the latter case survive much stronger rate of absorption, and for a small 
one they have a width mainly determined by the properties of the barrier itself. 

We also consider extension of our system to a periodic arrays of the sandwiches. We 
show that each of the resonances of the original system gives rise to a band, with the 
original resonance frequency at the center of the respective band. The band widths of the 
bands (and correspondingly their degrees of dispersion) are found to be determined by the 
frequency width of the respective parent resonances. This fact explains flattening photonic 
bands near ujt observed previously in numerical simulations of 2-D systems [14 ]. 
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work was supported by the NSF under grant No. DMR-9632789, by a CUNY collaborative 
grant, and by a PSC-CUNY research award. 
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CAPTIONS 

Figure 1. Transmission as a function of — . The ratio of — is equal to 3 for both plots, 
(a) Values of a and b are 2 and 8 respectively, (b) a = 2 and b = 20. The polariton gap lies 
between normalized frequencies 1 and 3. 

Figure 2. Transmission in the presence of absorption for a = 2 and 6 = 8. The dashed 
line corresponds to an attenuation coefficient ^ = 0.005, and the solid line depicts the 
transmission for -3- = 0.05. 

Figure 3. The bands arising in the polariton gap in the case of periodic arrangement of 
the considered system, (a) a = 2, b = 8, (b) a = 2, b = 20. The polariton gap is shown by 
dashed lines. Dispersion curves in the forbidden gap correspond to the resonant peaks from 
Fig. (1). 
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